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Motivated by rapid advances in plasma based accelerators, we propose a simple mechanism for 
the production of highly coUimated quasi monochromatic electrons beams. By considering Compton 
scattering of electrons (in the plasma) with virtual photons - through an effective interaction — we 
demonstrate angular collimation with a divergence less than 3 mrad, quasi monoenergetic nature 
of the beam and an energy gain of OilGeV / cm), in laser induced accelerators. 

PACS numbers: 41.75. Jv,52.38.Kd, 52.40.Mj,52.25.Mq 
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Collective effects in laser plasma interaction play a sig- 
nificant role in table top particle acceleration. Dawson 
and Tajima [l[ first showed, through computer simula- 
tions, that a high intensity laser shone on a dense plasma 
( no ~ lO^^cm"^ ) could accelerate the electrons in the 
plasma up to a GeV over a short distance of about a 
centimetre. Subsequently, significant progress has been 
made both experimentally and theoretically [2"10| ad- 
dressing, mainly, two issues: (i) acceleration to high en- 
ergies, and (ii) achieving a good collimation so that they 
can be of ut ility in high energy experiments and other 
applications [10|. While beam accelerations up to a few 
GeV had been achieved prior to 2004, the spatial diver- 
gence of the beam ranged from a few to tens of degrees 
with a Maxwellian energy spectrum [^-Q and a large 
energy spread (AE/E) ~ 100%. Monoenergetic beams 
with small spatial divergence were achieved first by three 
groups henceforth referred to as I, II and III respec- 
tively. In I and II collimation was achieved by carefully 
choosing laser parameters and plasma density whereas III 
used a preformed plasma channel to achieve the same. 
Unlike in the earlier experiments, the distribution was 
non Maxwellian in all the three cases. A few more 
salient points: In I the electron beam spectrum contained 
g - 2 ± 0.5 nC of charge with an energy E 170 ± 20 
MeV within the collimation angle A9 = 10 mrad. II 
reports Q ^ 22 pC of charge with A6 < 50 mrad. Fi- 
nally, III found Q - 0.32 nC of charge with - 86 ± 1.8 
McV and A9 ^ 3 mrad. Promising that these results arc 
- for future accelerator designs and experiments, theo- 
retical analyses have so far been computation intensive, 
employing PIC (three dimensional particlc-in-cell) sim- 
ulations [ll[ which have been used e.g., in [1, Q, and 
Vlasov models [l3|- Depending on the regime of interest, 
the acceleration mechanism is attributed to either laser 
wake field acceleration (in the linear regime), or stimu- 
lated Raman scattering (in the nonlinear regime 13[ . In 
this letter, in contrast to the above mentioned methods. 



we take an effective field theory approach to unravel the 
physics behind Plasma acceleration. We study accelera- 
tion via Compton scattering of photons by electrons in 
the medium, by explicitly incorporating the collective na- 
ture of the plasma at every level. We hope to show that 
the essential physics behind plasma accelerators can be 
captured by this approach and establish collimation and 
the monoenergetic nature of the scattered electrons. The 
attainment of high energies IGeV) over a distance 
of about a cm, will be addressed rather qualitatively, and 
will be treated more completely in a future work. 
The Effective Interaction: In the radiation gauge 
^ which we employ throughout, the constitutive 
Maxwell equations, which incorporate electric response 
of the plasma, have the form 



A'Kp{u!, k) ~ ~Loei{ui, k)k ■ A; 
[u:^ei{u:, k)Vl^ + {u:^et{u:, k) - P)rlAA, (1) 



where, Vl^ = hkj/k'^;Vlj = (6^, - Vl^) are the pro- 
jection operators for longitudinal and transverse modes 
respectively and e; and £t arc the corresponding per- 
mittivities. In the coUisionless limit and in the regime 
>> 1 relevant to us, we obtain the well known rela- 
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tions u}^ — a^t\k\ = ujp where, a; 
in terms of the plasma frequency LOp — AnriQe^ / m and 
the thermal velocity of the electron vt = ^jTjrn. T is 
the plasma temperature and m is the electronic mass. 

The lagrangian of the system including the medium re- 
sponses can be written as (we employ the units h = c = 

kB = 1), C = Cnirac + ^Int + ^Field, whcrC Coirac = 

Vj(a;)(i7^5^ - m)il){x); Lint = -e'ip{x)^'^'ip{x)A^{x) re- 
tain their usual forms while the field part gets modified 
to 



C 



where. 
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The corresponding hamiltonian for the field reads 



n 



Field 



kikj d{u!ei{uj, k) 



duj 



(3) 



The field energy can be interpreted as a sum of free har- 
monic oscillators provided that longitudinal and trans- 
verse components of the field are multiplied by respective 



normalization factors Nj' 



-^{u!'^et)\uj=ujj; IM- Using the dispersions relations we 



may further simplify the factors to Ni'' ~ ^t3"(1 — ^^'"^ 



2^ . Field quantization is then implemented via 
equal time commutation relations 



y), where U.i{r,t) = 
J d^r'dt'e^jit - t',r- r')diAj{t\r') The plane wave ex- 
pansion for the photons given by 
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where, A = 1,2 stand for transverse polarizations and A = 
3 denotes the longitudinal polarization. With an explicit 
employment of Ni^t leads to the standard commutation 



relations [ar, , at 



5{k - k')Sx,x> 



Once the quantization rules are established, it is 
straight forward to write Feynman rules (which we 
shall not do here), and also the amplitudes for any 
process. The relevant process here is Compton scat- 
tering j{uii,ki,ai) + e^{Ei,pi,Si) -f{ujf,kf,af) + 
e~{Ef,Pf,Sf) where a) and {E,p,s) refer to en- 

ergy, momentum and polarization for the photon and the 
electron respectively. Crucially, Compton scattering in 
the plasma medium involves longitudinal polarization of 
the photon in addition to the transverse modes. Accord- 
ingly, the scattering cross section is an incoherent sum of 
four contributions coming from the longitudinal (L) and 
transverse (T) modes in initial and final states. The spin 
averaged/summed scattering cross section for the process 
is given by da = ^ dai/; i, f = L, T. The target elec- 
trons belong to the plasma at a definite temperature, and 
this calls for a further averaging of the cross section wrt 
to the initial electron distribution. All the results pre- 
sented in this paper are reported after such an average. 
Thus the medium plays a dual role of (i) modifying the 
dispersion relations and the scattering amplitude through 
its permittivity and (ii) further contributing to the aver- 
age cross section. The role of the collective nature of the 
plasma is central to the analysis. It is pertinent to re- 
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FIG. 1: Angular regions of positive momentum transfer 
to the electron (in eV). The scattering angle 6e is 
defined wrt the beam direction. 



mark that the density and the temperature of the plasma 
have a highly nontrivial role to play. 

The basic quantities of interest are angular regions 
with momentum gain for the scattered electrons, the an- 
gular distribution {dN/d{cos6e)) and the corresponding 
energy spectrum (dN/dE); all of them are determined 
experimentally. The momentum gain can be determined 
from dispersion relations, and the other two quantities by 
multiplying the corresponding cross sections by the ini- 
tial photon flux. It is also of great interest to estimate yet 
another important quantity, the equivalent electric field 
that causes the acceleration, since it is a measure of gi- 
ant acceleration gradients which characterize all plasma 
based accelerators. We present the results below. 
Results and Discussion: The plasma parameters will 
be the same as chosen in the experiments. We take the 
plasma density uq = 6.5 x 10^^cm~^ and temperature 
To = 50 eV. The initial photon frequency uji ~ l.btJp, 
corresponding to a wavelength A — 9.31/im. The beam 
intensity is taken to be / = 3.2 x 10^^Wcm~-^ . First of 
all, consider the momentum transferred to the electron 
in each case as the scattering angle is varied: Fig [T] shows 
regions of positive momentum gain for each of the scat- 
tering processes. The angular width is very small, with 
A9 ~ 3 — 10 mrad with the TL and LL cases having 
the least divergence. We further note that the dominant 
contribution to AP comes from TL and LL scatterings 
which can be 10"^ larger than the other two cases. 
Collimation and Energy Spectrum: We now con- 
sider the number distribution {dN/d{cos9e)), in FiglJl 
for the electrons scattered in the angular regions with 
the positive momentum gain (See Fig. [T]); for clar- 
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FIG. 2: (a) Number density and (b) the corresponding 
cumulative distribution of scattered electrons as a 

function of Op 
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FIG. 4: (a)Plot of logj^Q Epff, equivalent electric field 
(V/cm) as a function of 9e- (b) The corresponding 
cumulative distribution. 
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FIG. 3: (a)Number density and (b) the cumulative 
distribution of the scattered electrons as a function of 
final electron energy Ef. Note the non Maxwellian 
nature. 




FIG. 5: (a) Extrapolated curve for dN/dEf after lO'' 
scatterings. (b)The corresponding cumulative 
distribution. 



ity, we show both the density as well as the cumula- 
tive distribution. Note again that the contribution to 
{dN/d{cos9e)) is dominated by LL scattering (FigH)). In- 
deed, the cross sections A.<7if - obtained by integration 
over the angular regions shown in FiglT] - are respec- 
tively given by AaLL = 1-2 x 10^, Acttl = 2.1 x 10^ 
AdLT = 1.3 X 10-^ Acttt = 0.0028 in units of = 
o? jrv? = 0.08 X 10"^^cm^, where tq is the classical 
electron radius. These numbers may be compared with 
the total cross section or ~ 10^° which is itself 10^ 
times larger than the Klien-Nishina cross section which 
is valid in free space. Fig 12] clearly demonstrates col- 
limation which is restricted to a small angular width 
8Q ^ 2 — 3mrad centered at Q^, = 1A74. This is in excel- 



lent agreement with the experimental result of Geddes 
et al [8| who report the same angular divergence, and 
is comparable to the findings of Faure et al [6| who re- 
port a divergence which is thrice the value found here. 
This establishes the first goal that we set out to demon- 
strate. It is important to note that the angle at which 
coUimation takes place depends on plasma parameters 
and the incoming beam energy. We report that we have 
verified that coUimation is stable against small perturba- 
tions in parameter values. Under large changes, e.g., in 
no, the beam quality is seen to deteriorate, as seen by 
experiments. The energy distribution of the electrons in 
the angular region around collimating cone is shown in 
FigUl The energy clearly peaks in the collimating cone, 
achieving a maximum density dN/dE = 2.5 x 10^"^ at 
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E = 19.46 cV with a small width AE/E = 2.37%. Note 
further that the distribution is non Maxwellian, as re- 
ported in experiments jfinSi]- 

Equivalent Electric Field: The equivalent electric 
field Eeff in this formalism is given by Eeff{9e) = 
= ^{AEl^), where nl is the mean pho- 
ton density with longitudinal polarization and El is the 
longitudinal component of the energy. Rather than esti- 
mate n'^ in terms of laser intensity and the medium prop- 
erties, we take the value of £'e//(^'e) as a phenomenologi- 
cal input from the experiments [l|, |6[ , who report a value 
w lO^Vcm"^ . This yields the value nj^ w 10^^ cm~^. 
At the collimating angle de = 1.47 rad, for the Eeff ~ 
0.8 X lO^Vcm^^ which compares excellently with the es- 
timate of [l| who find Eeff ^ lO^Vcm^^ . Employing 
this value, we plot Eeff{9e) in Fig|4]as a function of the 
scattering angle. The utility of this curve becomes clear 
below. 

Plasma Acceleration: Finally, we note that the energy 
of the scattered electron is peaked at £f ~ 20 eV, which 
is less than the thermal energy (50 eV). It might thus 
appear that the present efi^ective field theory approach 
fails to capture the essential feature of the plasma ac- 
celeration, viz., very large energies. We now argue that 
this is not necessarily a draw back: only a single scat- 
tering has been considered here while in reality the final 
energy owes its value to multiple scatterings. To put 
this on a more quantitative footing observe that the dis- 
tance traversed by the electron due to a single scattering, 
i ~ (nj^a)^^ w 0.3 nm where ri[ is taken from the above 
estimate. The number of scatterings in a distance of ^ 3 
mm is O(IO^), which would contribute to a cumulative 
energy £ ^ 100 McV. which is not far away from the ex- 
perimentally reported numbers. The result correspond- 
ing to this extrapolation is shown in FigE] for electrons 
accelerated over a distance of a centimeter. The energy 
spectrum shows several peaks and is remarkably akin to 
the electron spectrum unraveled in Q from deconvolu- 
tion of their experimental data. Note however that our 
data gives two closely spaced prominent peaks, one at 
£ Ri 118.5 MeV which corresponds at a scattering angle 
0e = 1.45 rad and the other at 5 « 150 MeV at 0e = 1.469 
rad. The second peak is relevant since it occurs in the 
collimating angular region. To conclude, we have seen 



that the physics behind collimation, quasi mono energetic 
spectrum and emergence of strong accelerating fields in 
laser plasma accelerators can be understood in terms of 
an effective lagrangian. This may be employed to study 
beat wave acceleration and plasma wake field accelera- 
tion with equal ease. A careful evaluation of multi scat- 
tering cross sections and an explicit incorporation of the 
beam properties would yield truly reliable quantitative 
estimates. Finally, the current approach allows a study 
higher order quantum corrections, similar to QED cor- 
rections to synchrotron radiation, which are not easily 
accessible in other methods. 
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